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METHOD FOR DETERMINING 
PRESSURE OF EARTH FORMATIONS 

FIELD OF THE INVENTION 
The invention relates to determination of properties of formations surrounding an 
earth borehole and, more particularly, to a method for determining properties including 
the leak-off rate of a mudcake, the perturbing effect of drilling fluid leak-off, and the 
undisturbed virgin formation pressure. 

BACKGROUND OF THE INVENTION 

A serious difficulty of formation pressure determination during drilling operations 
is related to the pressure build-up around a wellbore exposed to overbalanced pressure 
and subject to filtrate leak-off called supercharging. This pressure build-up is 
accompanied by filter cake deposition and growth externally, at the sand face, and 
internally due to the mud filtrate invasion. Thus, the filter cake hydraulic conductivity 
changes with time, affecting the pressure drop across it and therefore the pressure 
behind it, at the sand face. This makes it difficult to predict the evolution of the pressure 
profile with time, even if the history of local wellbore pressure variation has been 
recorded. 

Existing formation pressure measurements, made with so-called formation 
testing tools which probe the formations, often read high compared to the actual 
reservoir pressure far from the borehole, due to the supercharging effect. There are 
currently no known commercially viable techniques for the determination of the 
formation pressure in relatively low permeability reservoirs (below approximately 1 
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mD/cp) during drilling operations which adequately account for supercharging. The main 
difficulties are related to (1 ) the poor filter cake property, (2) the long actual time of 
wellbore exposure to overbalanced pressure, and (3) the practical time constraints, 
which require the pressure measurements to be carried out during a rather short time 
compared to the time of pressure build-up around a wellbore. These constraints make it 
difficult, if not impossible, to sense the far field formation pressure, at the boundary of 
the pressure build-up zone, with the usual transient pressure testing techniques, 
because of the slow pressure wave propagation inherent in low permeability formations. 

Accordingly, while existing tools and techniques can often work well in relatively 
high permeability formations, where supercharging easily dissipates, e.g. during tool 
setting, there is a need for a technique that can be successfully employed in relatively 
low permeability formations. It is further desirable to have a technique that is applicable 
to formations of wide ranging permeability, irrespective of the origin of the 
supercharging. There is also a need for accurate determination of filtrate leak-off 
parameters. It is among the objects of the present invention to address these needs. 

SUMMARY OF THE INVENTION 

In accordance with an embodiment of the invention, a method is set forth for 
determining the virgin formation pressure at a particular depth region of earth formations 
surrounding a borehole drilled using drilling mud, and on which a mudcake has formed, 
comprising the following steps: keeping track of the time since cessation of drilling at 
said depth region; deriving formation permeability at said depth region; causing wellbore 
pressure to vary periodically in time and determining, at said depth region, the periodic 
component and the non-periodic component of pressure measured in the formations 
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adjacent the mudcake; determining, using said time, said periodic component and said 
permeability, the formation pressure diffusivity and transmissibility and an estimate of 
the size of the pressure build-up zone around the wellbore at said depth region of the 
formations; determining, using said time, said formation pressure diffusivity and 
transmissibility, and said non-periodic component, the leak-off rate of the mudcake at 
said depth region; determining, using said leak-off rate, the pressure gradient in the 
formations adjacent the mudcake at said depth region; and extrapolating, using said 
pressure gradient and said size of the pressure build-up zone, to determine the virgin 
formation pressure. 

In accordance with a further embodiment of the invention, a method is set forth 
for determining the leak-off rate of a mudcake formed, at a particular depth region, on a 
borehole drilled in formations using drilling mud, and on which a mudcake has formed, 
comprising the following steps: deriving formation permeability at the depth region; 
causing wellbore pressure to vary periodically in time, and measuring, at the depth 
region, the time varying pressure in the borehole and the time varying pressure in the 
formations adjacent the mudcake; determining, at the depth region, an estimate of the 
flow resistance of the mudcake from the derived permeability and components of the 
measured pressure in the borehole and the measured pressure in the formations 
adjacent the mudcake; and determining, at the depth region, the leak-off rate of the 
mudcake from the estimated flow resistance and the measured pressure in the borehole 
and the measured pressure in the formations adjacent the mudcake. The virgin 
reservoir pressure can then be obtained by: determining, at the depth region, the 
pressure excess in the formations adjacent the mudcake from said derived permeability, 
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said leak-off rate, and said time since cessation of drilling; and determining, at said 
depth region, the virgin reservoir pressure from said measured pressure in the 
formations adjacent the mudcake and said pressure excess in the formations. 

Further features and advantages of the invention will become more readily 
apparent from the following detailed description when taken in conjunction with the 
accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a diagram, partially in block form, of a well logging apparatus that can 
be used in practicing embodiments of the invention. 

Figure 2 is a diagram of a downhole tool which can be used in practicing 
embodiments of the invention. 

Figure 3 is a diagram of logging-while-drilling apparatus that can be used in 
practicing embodiments of the invention. 

Figurer 4 is a graph of the quasi-steady pore pressure profile around the well 

bore. 

Figure 5 is a graph of dimensionless depth of pressure wave propagation into the 
reservoir. 

Figure 6 is a graph of formation response at the sand face. 

Figure 7 is a diagram of average pore pressure around a wellbore during pulse 
testing. Solid lines are shown, in the presence of pressure build-up; dashed lines, 
without build-up. 

Figure 8 is a graph showing pressure response at the wellbore to multiple-pulse 
production. 
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Figure 9 is a graph illustrating wellbore storage effect on pore pressure response 
at the wellbore for step-wise production for different ratios of formation to storage 
volume characteristic times. 

Figure 10 is a flow diagram of the steps of an embodiment of the invention. 

Figures 11 and 12 illustrate, respectively, testing in a pumping injection mode 
and in a production mode. 

Figure 13, which includes Figures 13A and 13B placed one below another, is a 
flow diagram of the steps of a further embodiment of the invention. 

Figure 14 shows graphs of the modulus (top track) and argument (bottom track) 
of the complex transfer function linking formation pressure at the sandface to wellbore 
pressure, plotted against frequency (in Hz). 

Figure 15 shows graphs of the modulus (top two tracks) and argument (bottom 
track) of the complex transfer function linking formation sandface pressure to wellbore 
pressure, as a function of dimensionless frequency co D = corl Ik, for a variety of values 

of mudcake skin. The upper two tracks repeat the same information, against linear and 
logarithmic y-axes. 

DETAILED DESCRIPTION 
Figure 1 illustrates a type of equipment that can be utilized in practicing 
embodiments of the invention. Figure 1 shows the borehole 32 that has been drilled in 
formations 31 , in known manner, with drilling equipment, and using drilling fluid or mud 
that has resulted in a mudcake represented at 35. For each depth region of interest, the 
time since cessation of drilling is kept track of, in known manner, for example by using a 
clock or other timing means, processor, and/or recorder. A formation tester apparatus 
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or device 100 is suspended in the borehole 32 on an armored multiconductor cable 33, 
the length of which substantially determines the depth of the device 1 00. Known depth 
gauge apparatus (not shown) is provided to measure cable displacement over a sheave 
wheel (not shown) and thus the depth of logging device 100 in the borehole 32. 
Circuitry 51 , shown at the surface although portions thereof may typically be downhole, 
represents control and communication circuitry for the investigating apparatus. Also 
shown at the surface are processor 50 and recorder 90. These may all generally be of 
known type, and include appropriate clock or other timing means. 

The logging device or tool 100 has an elongated body 105 which encloses the 
downhole portion of the device controls, chambers, measurement means, etc. 
Reference can be made, for example, U.S. Patents 3,934,468, and 4,860,581, which 
describe devices of suitable general type. One or more arms 1 23 can be mounted on 
pistons 125 which extend, e.g. under control from the surface, to set the tool. The 
logging device includes one or more probe modules that include a probe assembly 210 
having a probe that is outwardly displaced into contact with the borehole wall, piercing 
the mudcake 35 and communicating with the formations. The equipment and methods 
for taking individual hydrostatic pressure measurements and/or probe pressure 
measurements are well known in the art, and the logging device 100 is provided with 
these known capabilities. Referring to Figure 2, there is shown a portion of the well 
logging device 100 which can be used to practice a form of the invention wherein the 
variation in borehole pressure is implemented by the logging device itself (which, for 
purposes hereof includes any downhole equipment, wireline or otherwise) and is 
localized in the region where the device is positioned in the borehole at a given time. 
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(Reference can be made to U.S. Patent 5,789,669.) The device includes inflatable 
packers 431 and 432, which can be of a type that is known in the art, together with 
suitable activation means (not shown). When inflated, the packers 431 and 432 isolate 
the region 450 of the borehole, and the probe 446, shown with its own setting pistons 
447, operates from within the isolated region and communicates with the formations 
adjacent the mudcake. A pump-out module 475 which can be of a known type (see, for 
example, U.S. Patent No. 4,860,581), includes a pump and a valve, and the pump-out 
module 475 communicates via a line 478 with a borehole outside the isolated region 
450, and via a line 479, through the packer 431 , with the isolated region 450 of the 
borehole. The packers 431 , 432 and the pump-out module 475 can be controlled from 
the surface. The borehole pressure in the isolated region is measured by pressure 
gauge 492, and the probe pressure is measured by the pressure gauge 493. The 
borehole pressure outside the isolated region can be measured by pressure gauge 494. 
Embodiments hereof can utilize pumping and/or suction ports in the testing phase, and 
it will be understood that multiple pumping and/or suction ports can be provided. 

Embodiments of the present invention can also be practiced using measurement- 
while-drilling ("MWD") equipment (which includes measuring while tripping). Figure 3 
illustrates a drilling rig that includes a drill string 320, a drill bit 350, and MWD 
equipment 360 that can communicate with surface equipment (not shown) by known 
telemetry means. Preferably, the MWD equipment is provided with packers 361 and 
362. A device 365 is also shown, which includes probe(s) and measurement 
capabilities similar to the device described in conjunction with Figure 2. 
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The pressure build-up around the wellbore in relatively low permeability 
formations (such as k = 10" 1 mD) during drilling operations is a slow process, which 
usually lasts a few days and affects a relatively small neighborhood of the wellbore. The 
radius of the zone with elevated pressure around the wellbore can be estimated, using 
dimensional analysis. 

Assume that Darcy's law governs the flow in the reservoir 

v = --V P (1) 
M 

where v is the fluid flow velocity, // is the fluid viscosity and p is the pore pressure, 
which has to satisfy the pressure diffusivity equation 

dp __,2 fcS 

^r = 7 V 2 p, n = — (2) 
ot q>/x 

where t is the time, B is the bulk modulus of the rock saturated with fluid, ^ is the 
porosity and rj is the pressure diffusivity (see G.I Barenblatt, V.M. Entov and V.M. 
Ryzhik: Theory of Fluid Flows Through Natural Rocks, Dordrecht: Kluwer, 1990). 

If the time of exposure of the wellbore to overbalance pressure, t e , is known, 
then the radius of the zone with elevated pressure around it can be estimated as 

r e * 2^ (3) 
Using, for example, the following data: k = 10" 3 - 10" 1 mD , B = lGPa , // = lcp and 
<j> = 0.2 , one would obtain rj = (5 - 500) • 10" 6 m 2 /s. For the pressure build-up time 
t = 1 day, one finds 
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r e * 1.3 -13 m (4) 
The depth of investigation by conventional transient pressure testing, r f , also can be 
estimated, using the same formula (3). For example, if the investigation times are t f = 2 
hours, 20 min and 2 min, then the ratio rjr e can be respectively estimated as 

rjr e = y/We * 0.29, 0.12, 0.04 (5) 
This means that only first 29%, 12% and 4%, respectively, of the thickness of the 
pressure build-up zone can be sensed by the methods of transient pressure testing. 

The analysis of pressure build-up around the wellbore during drilling requires 
coupled consideration of the pressure wave propagation and the filter cake growth, 
induced by mud filtrate leak-off and usually restricted by the mud circulation inside the 
wellbore. If the overbalance pressure applied during drilling operations does not 
change dramatically, the transient pressure evolution around the wellbore can be 
approximated by the quasi-steady pressure behavior 



P{r,t) = 



Po> r > r e (t) 

where p 0 is the original formation pressure, p sf (t) is the pressure at the sand face, r w 
is the wellbore radius and r e (t) is the radius of zone around the wellbore with build-up 

pressure. The schematic of the pore pressure profile is shown in Figure 4. During the 
initial phase of wellbore exposure to overbalance, the pressure at the sand face, p sj , is 

equal to the wellbore pressure, p w . Then, the sand face pressure decreases with the 
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increase in the filter cake thickness and its hydraulic resistance due to the pressure 
drop across the filter cake, Ap = p w - p sj . 

If the filter cake permeability is small compared to that of the formation, the sand 
face pressure, p sf , falls quickly to the initial formation pressure, p 0 . If, however, the 

formation permeability is small and therefore the leak-off through the sand face is 
restricted, the filter cake is not built efficiently and the exposure of the formation to the 
overbalanced pressure can continue indefinitely. 

The unknown functions, p sf (t) and r e (t) , can be found from the pressure 

diffusivity equation (2) coupled with the model of the filter cake growth at the sand face. 
This analysis can be carried out for a simple model of the filter cake growth, based on 
the following assumptions: the porosity and permeability of filter cake are constant; the 
volumetric concentration of solids in the mud, filling the wellbore, is constant; the filtrate 
invading into the formation is fully miscible with the reservoir fluid; the filtrate viscosity is 
equal to that of reservoir fluid; and both spurt loss and internal filter cake formation are 
neglected. It is also assumed in this analysis that the filter cake permeability is much 
smaller than the reservoir permeability and the filter cake thickness, growing with time, 
is small compared to the wellbore radius. Under these assumptions, the flow through 
the filter cake can be considered as quasi-steady and one-dimensional at any time and 
therefore the pressure variation across the filter cake is linear as shown in Figure 4. 
The sand-face pressure, p sf (t) , is affected by a lot of factors, including the 

reservoir hydraulic conductivity, the leak-off rate and the rate of mud circulation. It also 
depends on the filter cake hydraulic resistance, which varies with time. Despite this 
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complexity, the boundary of the pressure perturbation zone, r e (t), plotted in appropriate 
dimensionless variables, is found to not practically depend on the filter cake growth 
dynamics and can be approximated by a universal function Z e (T), shown in Figure 5, 
where 



Z e (T) = Y(T)-l, y = p- 



T = £ (7) 



r w J 



Since the time of wellbore exposure to overbalance pressure, t e , is usually known, the 
only parameter, which is needed for the estimation of the radius of zone with perturbed 
pressure, r e (t e ) , is the pressure diffusivity, 77 , which is involved in the definition of the 
dimensionless time T . 

Assume that 77 has been found somehow and therefore the boundary r e (f e ) 



r e {t e ) = r\z e 



r 2 , 



+ 1 (8) 



Then, one has to measure the pore pressure at the sand face, p sf (t e ), and at an 

intermediate point r = r m inside the zone r w <r <r e (t e ) in order to find the formation 
pressure 

n _ P m jpgfcAj - Psf log(r e /r m ) 

p °- iSfcTy ' Pm=p(rJ (9) 

The sand-face pressure p sf (t e ) can be measured by currently available wireline testing 
tools and therefore, in order to obtain the formation pressure, p Q , one has to determine 
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the two parameters only - the pressure diffusivity, 77, and the pressure at some 
distance from the wellbore, p m , or alternatively the pressure gradient at the sand face 

vp s/ (t e )^ Pm(te) " Ps/(te) (10) 

Thus, if the formation transmissibility, kh/ju , which involves the interval thickness h , is 
known, the determination of the formation pressure, p 0 , is equivalent to the 
determination of the quasi-steady leak-off rate, q L (t e ) , at the end of the pressure build- 
up phase 

^0 = ^^vp s/ (t e ) (11) 

As shown below, q L , can be determined using pulse-harmonic tests, which can be 
carried out with appropriately chosen testing frequencies and pumping rates. 

In the following analysis of determination of far field formation pressure using 
pulse-harmonic testing, it is assumed that total testing time is small compared to the 
pressure build-up time (the time of borehole exposure to pressure overbalance); the 
pre-test volume is small compared to the total volume produced during testing, and the 
filter cake is tota ll y removed during pre-test. For simplicity, variation of the pressure 
diffusivity and the formation transmissibility versus the distance from the wellbore are 
ignored. 

Consider the situation just before pulse-harmonic testing, i.e. at t = t e . The 
pressure around the wellbore, p e (r) = p(r,t e ), specifies the initial condition with 
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respect to the testing time r = t - t e . Using the same notation for the pressure, p(r,r) , 
one has 



p(r,0) = p e (r), r>r w (12) 
As mentioned above, the function p e (r) is usually unknown except, its boundary value, 
Ptvo = P e ( r w) • which can be measured or estimated, using conventional formation 

testing. Using Eq. (6), the initial pressure profile around the wellbore before testing can 
be expressed as 

P.W = Po + ( P.-p 0 )^gg, r. SrS r.(0 (13) 

and the corresponding quasi-steady leak-off rate from the wellbore interval of thickness 
h is 

_ 2xkh p w0 -p 0 
q ^ M log[r e (t e )/rJ (M) 

This leak-off rate, q L , is unknown in advance and its determination would be equivalent 
to the determination of the two parameters: the radius of the pressure build-up zone, 
r e (t e ) , and the formation pressure, p 0 . 

Using Eq. (14), the initial pressure profile can be represented in the equivalent 

form 



P e (r) = Puo-Pz.log — 



V 



^ = ^wT (15) 
2nkh 
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Generally speaking, the parameter q> L could be determined, using, for example, the 
conventional pressure build-up technique, if one could seal instantaneously the sand 
face of the wellbore interval and monitor the pressure relaxation, p w (r) , behind the 

sand face with time. Indeed, due to the superposition principle, the pressure response 
at the sealed sand face to the step-wise variation of the flow rate can be expressed as 

¥ w (t) = p w (t) - p w0 = -<p L F 0 (ijr/rZ) (16) 
Here, the function F 0 (a), where a = rjr/r^ , is given by the well-known solution of the 

pressure diffusivity equation (see, for example, H.S. Carslaw and J.C. Jaeger: 
Conduction of Heat in Solids, 2 nd Edition, Oxford: Clarendon Press, 1959) 



7T * 



4a? (l-e^)dg 



(17) 



where the J, and Y, are Bessel functions of the first and second kind, respectively, of 
order /, /' =0, 1 , and it is shown in Figure 6, reproduced from Carslaw et al., supra. 
Since, at large time 

f 2.2577/ 



¥ w (r) * -<p L log 



r 2 



(18) 



one could determine the two parameters, <p L and 7/r 2 , by plotting y/ w {r) versus logr. 

This straightforward approach, which is widely used in the well testing technology 
(see T.D. Streltsova: Well Testing in Heterogeneous Formations, Exxon Monograph, 
John Wiley and Sons, 1988), is, however, rather difficult to implement in reality. There 
are a few reasons for this. First of all, the necessary testing time in low permeability 
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formations is usually extensive . Secondly, the initial leak-off rate in a low permeability . 
formation is typically very small and can be very difficult to measure. The sealing of the 
sandface and the pressure monitoring is preferably done with great care so as not to 
disturb the formation and the pressure at the sandface. It is worth noting also that the 
sealing of the wellbore surface could be replaced by the pressure relaxation procedure, 
which would prevent the leak-off, but this is not much easier to implement because the 
detection of a very small leak-off can be even more challenging. Thus, a different type 
of pressure testing procedure is needed. Pulse-harmonic testing has the advantage of 
not compromising the accuracy of measurements and the amount of information to be 
extracted from the data is comparable to that, which maybe extracted by conventional 
means. 

Consider the pressure evolution around the wellbore during pulse-harmonic 
testing with a production rate q w (r) , having a period f. Using the superposition 
principle, one can represent the production rate perturbation during testing, 
<?(*") = <lw ( T ) + Ql , as a sum of its periodic component with zero average rate, q p (r) , 
and the constant average rate, q a , i.e. 

<?(*") = <Z P ( T ) + 9a. q a =q w + q L , q p (r) = q w (r)-q w (19) 

where 




The unknown leak-off rate, q L , has been added to the production rate q w (r) to 
compensate for the initial non-uniform pressure profile (15) around the wellbore. The 
advantage of this testing procedure is that the periodic part, q p (r) , can be tuned for 
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different depths of investigation, R « 2^nrfT , by changing the angular frequency 

o) = (see Stretsolva, supra). The testing time is comparable with the period f 
and is usually much shorter than the duration of a pressure build-up after shut-in. At the 
same time, the average rate, q w , should not depend too much on the characteristics of 
the hardware (pumps, pressure gauges, flow meters). It can be tuned by choosing, for 
example, appropriate amplitudes, q 0 , and durations, t 0J of production pulses and the 

ratio t 0 /f (Figure 8). The interpretation of the responses to the periodic component, 
q p (r) , and non-periodic component, q a , of the production rate then can be carried out 
independently. 

The other advantage of this superposition is that the periodic component, q p (r) , 
does not involve the unknown initial leak-off rate, q h% and the extraction of the pressure 
response to the periodic rate q p (r) , from the measured pressure variation at the 
wellbore, y/ w (r) , is a standard task in the practice of pulse-harmonic testing (see 

Streltsova, supra). Processing the pressure response to the periodic component, allows 
one to determine the pressure diffusivity, 77, and the formation transmissibility, kh/ju . 
This reduces the number of unknown parameters in the presentation of the initial 
pressure profile before testing, determined by Eqs. (13) and (8), to only one - the 
formation pressure, p 0 . 

The determination of p 0 requires the processing of the wellbore pressure 
response to the non-periodic component of the production rate, which is represented by 
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the average constant rate, q a . Using the superposition principle, this response can be 
expressed similarly to (16) as 



r.M-«. +ft )MW*). < 21 > 

Here, ^ a (r) is the measured pressure response minus the periodic component; the 
parameter p w is already known, and the parameter <p L is still unknown. 

The function F 0 (a) is defined by (17) and shown in Figure 6. Since the 
pressure diffusivity, ?j, has already been determined from the pressure response to the 
periodic component, the argument a = rjr/r^ can be calculated. Now, compare Eq. (16) 
and Eq. (21 ). Eq. (16), which corresponds to the standard pressure build-up test, 
involves two unknowns, y L and rj, whereas Eq. (21) involves only a single unknown 
parameter, <p L . This advantage can be exploited to full extent. Indeed, the parameter <p L 
can be estimated, using the pulse-harmonic testing data, as 

<Pl=~V w - J, a{j }^ (22) 

Thus, the last term in the right-hand side of Eq. (22), which formally depends on the 
testing time r , has actually to be constant. This term can be estimated, using the 
pressure measurements in the wellbore, y/ a (r) , and the function F 0 (a) , representing 
the dimensionless reservoir pressure response to an average step-wise production rate. 

After the determination of the parameter <p L , the desired formation pressure can 
be estimated as 
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Po = Pu,o ~ <Pl log[r.(t,)/rJ (23) 
Eq. (22) can be also interpreted as follows. In the absence of the initial pressure build- 
up and the corresponding leak-off rate, the last term in its right-hand side has to be 
equal exactly to p w . This means that the difference between the two terms at q L * 0 

represents the effect of the "boundary condition" at the virtual moving boundary, 
corresponding to the pressure wave, propagating into the formation, as shown in Figure 
7. Here, the pressure profiles are plotted in the logarithmic scale I = logr for three 

sequential testing times r, < r 2 < r 3 . Since the average production rate is constant, the 
solid lines, representing the pressure profiles in presence of the initial pressure build-up, 
Pwo ~ Po > nave tne same slopes. The dashed lines represent the pressure profiles, 
which should be observed in the absence of the initial pressure build-up. It is assumed 
also that the velocity of the virtual front of the pressure wave, l = l M , propagating into 
the formation, is not affected by the pressure build-up. For this reason, the difference 
between the wellbore pressure behavior in the two cases is accumulated with time: 
Ap, < Ap 2 < Ap 3 . This accumulated difference makes the term -y/ a (r) = p wQ - p w (r) , 

involved in Eq. (22), larger than the denominator F 0 (rjr/r^j , which represents the 
response to the step-wise rate, y w , corresponding to the uniform initial pressure profile. 

In the following example, consider the multiple-pulse testing procedure, 
illustrated in Figure 8, with the production pulse amplitude q 0 , the production pulse 

duration t 0 , the period f and the time lag between two sequential pulses t, = f - 1 0 . 
The average production rate, q w , can be found from (20) as 
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q w = q 0 (to/T) (24) 

Using the superposition principle, the pressure response to the first production pulse at 
the wellbore can be represented as 



¥ w to = -? w |>o («) ~0(r- t 0 )F 0 (a,)] 
where 0(t) is the Heaviside unit step function and 



(25) 



Vw 2xkh 



ft \ 

If j 



r r 



(26) 



Using the measurements of the pressure perturbation at the first shut-in (the point A in 
Figure 8) and at the beginning of the second production period (the point B), y/ A and 

y/ B , the equation for the pressure diffusivity tj can be obtained 



Vb F 0 (fjT/r^)-F 0 (TjtJr^) 
After 7] has been found, the formation transmissibility can be calculated as 



(27) 



kh 


_ <7o 


ft > 






H 


2ny/ A 


If J 







(28) 

Now, the pressure response at the wellbore to the non-periodic rate, y/ a (r) , has to be 

extracted from the measured pressure curve OABCD... as shown in Figure 8. This 
means that at least the first three production pulses preferably should be involved in 
interpretation to allow the determination of y/ a (r) with confidence. Finally, the 
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parameter y L , which is proportional to the initial leak-off rate q L , can be found, using 
Eq. (22), and then the formation pressure is calculated from Eq. (23) 



Po=Pwo-PlI°& 



r w 



k 




+ 1 




r 2 





(29) 



where the function Z e (T) is shown in Figure 5. 

The graphical interpretation in Figure 7 aids in the understanding of the 
requirements of the pulse testing design, which should reduce possible interpretation 

errors. It is obvious that the average production rate q 0 (t 0 /f) should not be too high 

compared to the leak-off rate, otherwise the right-hand side of Eq. (22) will be small 
compared to the terms involved in their residual and therefore errors of their 
measurements may affect the accuracy of calculation of <p L . The best resolution should 

be achieved when q 0 (tjf) is close to the leak-off rate. In this case, the slopes of the 

local transient pressure profiles and the build-up pressure profile are equal but have 
opposite signs. 

The fluid volume, located between the pump and the wellbore surface (or sand 
face), which is known also as a storage volume, can distort the production pulses 
created at the pump. As a result of this distortion, the boundary condition at the 
wellbore surface does not match exactly the production schedule, generated by the 
pump, and therefore the pressure response is different from the obtained solution. This 
phenomenon, known as a wellbore (or tool) storage effect, can be important if the 
storage volume is large compared to the total production volume per testing cycle. 
Indeed, the storage volume is decompressed during production and pressurized during 
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injection cycles, damping the rate variation, induced by the pump, and therefore 
smoothing the formation response to it. If the compressibility of the fluid in the storage 
volume is constant, the storage effect can be investigated, using the Laplace 
transformation technique (see Barenblatt et al., supra, and Carslaw et al., supra). 

The fundamental solution for the step-wise production rate with amplitude q 0 and 
zero initial conditions is given (Carslaw et al., supra) by the formulae 

<p(r,T) = -<p 0 F s (a), a = ^, <Po = £^ (30) 
4a? ll-e-f)d4 

u(z) = yzJ 0 {z)-J 1 (z), v(z) = r zY 0 (z)-Y l (z) (32) 
It involves the additional dimensionless parameter y , which is determined as 

H =f%L, r P =& (33) 

which is the ratio of the two characteristic times, t s and r F , corresponding to the 
storage volume and the formation respectively. Here, is the storage volume and c 0 
is the fluid compressibility, which correlates the variation of the storage volume, AVg, 
with the pressure variation, Ap, as AV S = -c 0 V s AP. The solution (31)-(32) becomes 

identical to (17) at y = 0 . The function (2^) _1 F s (a) versus log 10 (a) for y' 1 = 0.5 , 1 , 

2, 4 and » is shown in Figure 9 (reproduced from Carslaw et al.). One can see that the 
storage effect is more pronounced at small time, especially for large y. This solution 
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can be used for the interpretation of the pulse testing data as outlined above instead of 
the solution (16)-(17). 

It will be understood that the described technique can be expanded to take into 
account the variation of the formation properties, i.e. the pressure diffusivity and 
transmissibility, with the distance from the wellbore due to invasion of mud filtrate into 
the formation during drilling. Pulse-harmonic testing with different frequencies can be 
used to discriminate the responses of the damaged zone and the undamaged 
formation. The design of the testing procedure in such a case would require some a 
priori information (at least, an order of magnitude estimate) about the formation 
transmissibility and diffusivity. If they vary significantly with distance from the wellbore, 
the interpretation of the pressure response to a non-periodic component of the 
production rate would need to be modified, and a longer testing time would generally be 
necessary. 

Figure 10 is a flow diagram of steps for practicing an embodiment of the 
invention, as described. The block 1003 represents keeping track of the time since 
cessation of drilling at the depth region(s) of interest. A pretest is performed (block 
1005) and downhole parameters, including permeability, are measured in conventional 
fashion (block 1010). Borehole pressure in the zone is increased (block 1020), and 
oscillated flow rate (block 1030). As discussed, the pressure can be controlled, for 
example, from the wellhead or between the dual packers. A first set of downhole 
parameters is determined (block 1040). In the present embodiment, this includes 
determining, using the periodic component of the measured pressure, the formation 
pressure diffusivity and transmissibility, and an estimate of the size of the pressure 



22 



build-up zone around the wellbore. Then, as described, this set of downhole 
parameters, and the non-periodic component of the measured pressure, are used to 
determine the filtrate leak-off rate and/or the pressure gradient (block 1060). The 
formation pressure can then be determined by extrapolation (block 1075). 

Figures 11 and 12 illustrate testing in a pumping/injection mode (Fig. 11) and a 
production mode (Fig. 12). 

For the pumping/injection mode of Fig. 1 1, a primary purpose is measurement of 
the hydraulic conductivity of the mudcake, which should not be significantly damaged, 
removed or modified if fluid is pumped through it into the formation. The packed off 
interval may be used to: a) reduce the effects of tool storage, b) selectively isolate a 
specific depth region for testing and/or c) to increase the surface area and to maintain 
an appropriate injection rate that will induce measurable pressure response behind the 
mudcake without formation fracturing, among others. In Fig. 1 1 , the time scale starts 
from the tool setting and probe penetration through the mudcake followed by the small 
volume pretest (shown at (a)) in order to cleanup the probe-formation interface and to 
establish good hydraulic communication between the pressure gauge (e.g. 493 in Fig. 
2) and the formation sand face. After pressure build-up (shown at (b)), the fluid is 
injected into the formation through the packed off interval covered by mudcake using 
pulses (shown at (c)), creating transient pressure response behind the mudcake. The 
pressure at the sand face measured with the probe increases during injection pulses 
and relaxes between them, whereas the interval pressure is maintained constant during 
injections. The two pressures measured by gauges 492 (interval) and 493 (probe) allow 
for the calculation of the mudcake hydraulic conductivity, as described below. It is 
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possible, using known methods to determine the diffusivity and the storativity 
respectively by employing low frequency and relatively high frequencies. 

For testing in a production mode, as illustrated in Fig. 12, the purposes include: 
(1 ) determining formation parameters (the pressure diffusivity and the pressure 
transmissibility or kh/|i ) using the periodic pressure response at the sand face to 
production pulses, and then (2) estimating the initial leak-off rate from the wellbore into 
the formation using the non-periodic pressure response. The analysis has been set 
forth in detail above. As shown in the Figure 12, the pre-test (a) is performed for 
mudcake cleanup and establishing good hydraulic communication between the tool and 
formation, followed by a few production pulses. The number of production pulses is 
preferably at least three. More pulses will tend to increase the resolution of the non- 
periodic part of the pressure response. 

A further embodiment of the invention will next be described, this embodiment 
including a technique for estimating the parameters of the mudcake which control filtrate 
leak-off rate, and for using this estimate in turn to estimate the true reservoir pressure 
from the measured sandface value. A flow diagram of the steps for practicing this 
embodiment is shown in Figure 13. 

The time post-drilling is kept track of (block 1 103). As represented by block 
1 105, a formation pressure measurement tool is deployed in the well, and set on the 
formation of interest. An estimate of the formation permeability is made (block 
1110). This can be done using standard means; for example, interpretation of pre- 
test pressure transients. This is combined with an estimate of the formation total 
compressibility, to obtain an estimate of the formation pressure diffusivity (block 
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1115). The wellbore pressure is caused to vary periodically in time (block 1 125) with 
significant frequency content in an appropriate frequency range, as discussed 
above, and treated further below. The time-varying pressures measured by the 
formation probe pressure sensor, and a pressure sensor in the wellbore (Figure 2), 
are measured and recorded (block 1 130). The time-periodic parts of the wellbore 
and formation pressure measurements are analyzed, using also the information on 
the formation permeability obtained from the pre-test, so as to give an estimate of 
the flow resistance of the mudcake (block 1 140). 

The estimated flow resistance of the mudcake is then combined with the measured 
wellbore and sandface pressures to estimate the filtrate leak-off rate (block 1 150). 
Then, as represented by the block 1160, the filtrate leak-off rate is combined with the 
estimated formation permeability and the time of exposure of the formation post- 
drilling, to estimate the pressure excess at the sandface due to leak-off (i.e. 
supercharging). This pressure excess is subtracted from the measured pressure, to 
yield an estimate of the true reservoir pressure uncontaminated by supercharging 
(block 1170). 

Further detail of the routine for this embodiment will next be described. 
Regarding step 1125, once the tool's probe is set and in pressure communication with 
the formation, steps are taken to induce modest amplitude, time periodic, absolute 
pressure variations within the wellbore, so as to create (a) measurable pressure 
disturbance within the wellbore at the tool, and (b) a measurable response to this 
disturbance, as seen by the pressure sensor in communication with the formation 
through the probe (e.g. Figure 2). 



25 



The wellbore pressure can be written as p w {t) = p w + w(p w (a>)e""), where p w 
denotes the (constant) background wellbore pressure about which the fluctuations take 
place, W() indicate the "real part" of argument, p w denotes the amplitude of the 

oscillation, a is the frequency. Mechanisms for generation of pressure variations within 
the formation include the response to changing filtrate loss rates through the mudcake 
(although other processes could contribute, e.g. elastic deformations of the rock or 
deformation of the mudcake itself). The frequency of the wellbore pressure fluctuations 
should be chosen so that the measured attenuation of pressure fluctuations across the 
mudcake is adequately sensitive to the flow resistance of the mudcake. Computed 
pressure responses are shown in Figures 14 and 15, and inspection of these indicates 

that a good choice of frequency is in the range a> D = a>r* /rj = o(l0" 2 to 10°) , because 
responses are not too small, nor dimensional frequencies too low [r w is the wellbore 
radius measured on the rock side of the mudcake, rj is the diffusivity of pressure within 
the formation, ando is the angular frequency of the induce pressure pulsations). 
Selection of frequency was treated above. A further consideration in selection of 
frequency is that it should be low enough that the depth of penetration of pressure 
disturbances is greater than the thickness of the mudcake, and this translates into the 
requirement that (f> c pc c cod 2 lk c «\, where d is the mudcake thickness, c c is the 
mudcake compressibility, <j> c is the mudcake porosity, fc c is the mudcake permeability 
and k c /<p c juc c is a measure of the diffusivity of pressure within the mudcake. 
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Regarding interpretation of attenuation of pressure fluctuations for the mudcake 
skin, the complex amplitude of axisymmetrical time harmonic pressure fluctuations 
within the formation, having angular frequency co, satisfies 



. « 1 d 

icop 

r dr 



dr j 



(34) 



where actual pressures are given by p(r,t) = W.(p{r,co)e i0 *), rj = k/<f>^c t , where k is the 
formation permeability, </> the formation porosity, ju the viscosity of the fluid in the pore 
space and c, the compressibility of the fluid-solid system (formation saturated with 
fluid). Pressure fluctuations decay at great distances, so p(r,co) -> 0 as r -» oo . At the 
wellbore wall, the mudcake is modeled as an infinitesimally thin "skin", across which 
there is a pressure loss proportional to the instantaneous flow rate, so that 



- P(r w , <o) = -rJS ^ (r w , co) , 
dr 



(35) 



where the non-dimensional parameter S is the standard skin factor familiar in well 
testing. It can be shown that 



P(r w ,a>) = p w {co)- 



fico 



+ 



f rr- \ 
ico 



(36) 



where the K's are modified Bessel functions, and the branch of the square root is 
chosen so as to ensure decay of pressure perturbations at large distances. 

Figures 14 and 15 show graphs of the modulus and argument of p{r w ,a>)/p w (co), 

as given by the above formula, plotted versus co or co D = ar^ / tj for a variety of values 
of S . In Figure 14, the formation permeability is 10mD, the porosity 20% of the 
formation fluid viscosity ImPa.s, the total compressibility 10" 8 Pa* 1 , the wellbore radius 
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0.1m, and the mudcake skin S=99.49 (corresponding to a cake of thickness 1mm and 
permeability 0.001 mD). For such a mudcake, the fluid loss rate driven by a 1 0Opsi 
pressure differential is 6.8x1 0" 5 cm/s. From Figure 1 5, it can be seen that if the values 
of 7 , o) and r w , and hence cd d , are known, then it is possible to estimate the value of 
S from the measured value of the ratio of the amplitudes of the sandface and wellbore 
pressure fluctuations, \p{r w ,co)/p w (©)| . In the present embodiment, the values of p w (a>) 

and p{r w ,co) are obtained from the measured time series of p w (t) and p(r w ,t) using 

standard signal processing methods. 

As a further refinement, the drilling fluid circulation rate and/or long-time average 
wellbore pressure can also be varied. Changes in circulation rate will cause erosion (or 
further growth) of the mudcake, and changes in filtration pressure will cause the cake to 
compact (or expand slightly). The cake skin at each circulation rate or overpressure 
can be estimated using the method just outlined, and by this means a table of values of 
S versus circulation rate (denoted as y) and/or filtration pressure (p w -p(r w ,t), 
denoted as Ap) can be created. The values stored in this table can be used in the step 
of block 1 1 50 (treated further below), so that the value of S corresponding to the 
current circulation conditions is used when evaluating the leak-off rate. Interpolation 
between measured values may be used. 

Regarding the step of block 1 150, the instantaneous pressure drop across the 
mudcake is related to the sandface pressure gradient by 

p w {t)-p{r w ,t) = -rj(r{t\ty(t))^ (r„,l). (37) 

and using Darcy's law at the sandface, 
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--%.<>1. (38) 
fx dr 

to relate the sandface pressure gradient to the filtrate leak-off flux, q , one obtains 

Using this expression, under the assumptions that (a) the fluid loss can be adequately 
described by the skin parameter S estimated above, and (b) sufficient data has been 
collected in the previous steps to permit extrapolation and interpolation to estimate S 
over the range of wellbore flow rates and pressures occurring between first exposure of 
the formation and the formation pressure measurement (or have a mechanistic model to 
link values of S measured at one set of wellbore conditions to those pertaining at 
another), the filtrate loss rate q(t) can be estimated given the measured time histories of 
wellbore and sandface pressures, p w (t) and p{r w ,i), respectively and information on 
the drilling fluid circulation rate. 

Regarding steps 1160 and 1170, the sandface pressure is related to the fluid 
leak-off rate through the familiar convolution integral 

^.0 = A.+JG(<-'W')* , . (40) 

'o 

where t 0 denotes the time at which the formation was first drilled, p 9 is the reservoir 
pressure at great distances from the well, G is the formation impulse response which 
contains as parameters the formation permeability {k) and pressure diffusivity (77), and 
q(f) is the filtrate leak-off rate time history estimated as described above. The 
functional form of G is well known in the art. 
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By comparing the predicted sandface pressure, given by the previous equation, 
with the sandface pressures actually measured, can be estimated. Stated another 



way, the quantity \G(t-t')q(t')dt' can be taken as an estimate of the overpressure due 



to supercharging, and subtracted from measured pressures so as to give an estimate of 
the true formation pressure. It will be understood that this embodiment relies on an 
indirect estimation of overpressures from filtercake resistance which affects the 
accuracy of the technique. The interpretation model assumes that that mudcake is thin, 
and behaves like a simple additional resistance to fluid flow between wellbore and 
formation. The technique may be modified to take account of the finite thickness of the 
cake, unsteady pressure diffusion within the cake itself, and/or interactions between the 
hydraulic properties of the cake and the changing wellbore pressure. 

While the invention has been described with respect to a limited number of 
embodiments, those skilled in the art, having benefit of this disclosure, will appreciate 
that other embodiments can be devised which do not depart from the scope of the 
invention as disclosed herein. For example, embodiments of the invention may be 
easily adapted and used to perform specific formation sampling or testing operations 
without departing from the spirit of the invention. Accordingly, the scope of the invention 
should be limited only by the attached claims. 
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